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At the LHC, the processes with several particles in the final state occur routinely. These processes 
provide a new domain to test the Standard Model of particle physics. Due to large gluon lumi- 
nosity at the LHC, gluon initiated processes will be quite important. We study the production 
of multi-vector bosons via gluon fusion. In particular we consider the process gg — > gyZ and 
compute its cross-section at hadron colliders. 
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1. Introduction 

The Large Hadron Collider (LHC) has begun its operation. The Standard Model (SM) seems 
to be in excellent agreement with the collected data. There has been searches for the hints of 
physics beyond the SM such as supersymmetry, large extra dimensions, etc. But, as of now, there 
is no clear evidence [|J]. The search for the Higgs boson is also continuing. With the lack of signals 
for beyond the SM scenarios, there is a need to look for the SM processes that were not accessible 
at the Tevatron. Most of such processes have several particles in the final state, and/or occur at the 
higher order. Such processes can also contribute to the background to the new physics signals. One 
such class of processes is multi vector boson production in association with one or more jets. 

At the LHC centre of mass energy, the collider has another useful feature. In the proton- 
proton collisions, the gluon luminosity can be quite significant. It can even dominate over the 
quark luminosity in certain kinematic domains. Therefore at the LHC, loop mediated gluon fusion 
processes can be important. Di-vector boson production via gluon fusion have been studied by 
many authors [pL 51 0L BL EJ. We consider a class of processes gg — > VV'g, where V and V' can 
be any allowed combination of electroweak vector bosons. At the leading order, these processes 
receive contribution from quark loop diagrams. The prototype diagrams are displayed in Fig ([!]). 
The computation for the process gg — >■ yyg has already been performed [0, |p. In this contribution 
to the proceedings, we particularly focus on a specific process gg — >• gyZ. We report our results of 
the cross-section calculation for this process at hadron colliders. After this brief introduction, in 
the next section, we describe our process and discuss the structure of the amplitude. In the section 
3, we provide details on our calculation and numerical checks. In the section 4, we present brief 
numerical results. We conclude in the last section. 




Figure 1: The prototype diagrams for the processes gg — > VV'g. The wavy lines represent the appropriate 
combination of the y, W, or Z boson. The last two classes (c) and (d) are relevant to WWg production only. 
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2. The Process 

The process gg — > gyZ contributes to the cross-section of the process pp — > yZ+ jet. This 
process can be a background to new physics scenarios, in particular to the technicolor scenarios. 
Although the process gg — > gyZ occurs at the one loop level, its cross-section can be significant 
due to large gluon luminosity at the LHC. This one-loop process receives contribution from two 
main classes of quark loop diagrams - pentagon and box type, as shown in (a) and (b) of Fig 
Box class of diagrams are due to the triple gluon vertices and can be further divided into three 
sub-classes. This sub-classification has it's own physical importance as we will see later. Other 
diagrams can be obtained by suitable permutation of external legs. For each quark flavour, there 
are 24 pentagon-type and 3 x 6 (= 18) box-type diagrams. Due to Furry 's theorem, only half of 
the 42 diagrams are independent. We work with all six quark flavours. Except ?-quark we treat all 
other quarks massless. Our one loop process, being at leading order, is expected to be finite. 

The amplitude for our process has very nice structure. Let us first consider the prototype box 
diagram. Since Z-boson has both vector and axial-vector couplings with quark, one can write this 
diagram as a sum of two pieces - vector piece and axial- vector piece 

Here i is the diagram index. It takes value from 1 to 18, for each quark flavour. Furthermore, q is 
the flavor index for the quark in the loop and runs from 1 to 6. 

The axial-vector piece will cancel when we sum over all box diagrams. This happens because 
we can categorize 18 box diagrams for each quark flavour in two parts. One part has momentum 
flowing in clockwise direction and other part has momentum flowing in the anticlockwise direction. 
These two parts have identical vector pieces, but the axial vector pieces have opposite signs. It leads 
to the cancellation of axial vector pieces. Therefore, the axial-vector coupling of the Z-boson to 
the quarks does not make any contribution to the full box amplitude. There is pairwise cancellation 
between charge-conjugated fermion loop diagrams as shown in Fig @. It means 

18 

E<s=0- (2-2) 
i=i 

The situation is different for the pentagon diagrams. We can again write the contribution of 
each pentagon diagram as the sum of vector and axial-vector pieces 

<p = <p+<p- (2-3) 

Here i is the diagram index. It runs from 1 to 24 for each quark flavour. The flavour index q runs 
from 1 to 6. Without color factor, the £f*'p for the diagrams with opposite direction of the flow of 
momentum have opposite signs, while the have same signs. Because of the opposite signs, 
the sum of the vector pieces of all pentagrams have same color factor, as the sum of box diagrams. 
This color factor is antisymmetric (i.e., it is proportional to the structure constants of the SU(3) Lie 
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Figure 2: Implication of Furry 's theorem for charge-conjugated box-diagrams. 

algebra, f a bc)- The color factor for the axial-vector piece is symmetric (i.e. proportional to d a bc)- 
This has been shown in Fig (0). Thus the full amplitude has following general structure 
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Figure 3: Implication of Furry's theorem for charge-conjugated pentagon-diagrams. 
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Mv^ are color ordered amplitudes for the vector and axial- vector parts of the amplitude, s/^ and ^ 
include full box and pentagon contributions for all six quark flavours. My\ and My\ represent light 
quark and ?-quark contributions respectively. The structure of the amplitude suggests that the vec- 
tor and axial-vector contributions should be separately gauge invariant. Moreover due to the color 
structure when we square the amplitude, the interference between the axial and vector contribution 



V,A 
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vanishes, i.e. 



M(ZZ^zyZ)\ 2 = [6\M V \ 2 +™\M A \ 2 1 . (2.8) 



Therefore the cross-section of the process is an incoherent sum of the vector and axial-vector 
contributions. Hereafter we will consider only the vector contribution, My. 

3. Calculation and Numerical Checks 

For each class of diagrams, we write down the prototype amplitudes using SM Feynman rules. 
The amplitude of all other diagrams are generated by appropriately permuting the external momenta 
and polarizations in our code. The quark loop traces are calculated in FORM in ^-dimensions [0]. 
The amplitude contains tensor loop integrals. In the case of pentagon-type diagrams, the most 
complicated integral is rank-5 tensor integral (E^ vpaS ); while for the box-type diagrams, rank-4 
tensor integral (D^ vpa ) is the most complicated one 

gtivpoS _ f d n k k^k v k p k a k 5 
J (27t) n NoN 1 N 2 N 3 N 4 ' 

J (2k)" N N l N 2 N 3 ' 

Here, Nj = kf — m 2 + ie and is the momentum of the i th internal line in the corresponding scalar 
integrals, n = (4 — 2e) and m q is the mass of the quark in the loop. Five point tensor and scalar 
integrals are written in terms of box tensor and scalar integrals using 4-dimensional Schouten Iden- 
tity. For example, a five point scalar integral, in 4-dimensions can be expressed in terms of five box 



scalar integrals [11] 



£b(0, 1,2,3,4) =£ciD£°, (3.3) 

i=0 

where D$ is the box scalar integral obtained after removing the i th propagator in Eq. It can be 
argued that any 0(e) correction to the above relation in n = (4 — 2e)- dimensions does not con- 
tribute to our amplitude. The box tensor integrals are reduced into the standard scalar integrals - 
Aq, Bq, Cq and Do using fortran routines ^ that follows from the reduction scheme developed by 



Oldenborgh and Vermaseren [10]. Thus we require box scalar integrals with two massive external 



legs, at the most. The scalar integrals with massive internal lines are calculated using OneLoop 



library Q12[]. Because of very large and complicated expression of the amplitude, we calculate the 
amplitude numerically before squaring it. This requires numerical evaluation of the polarization 
vectors of gauge bosons. We choose real basis, instead of helicity basis for the polarization vectors 
to calculate the amplitude. This is to reduce the size of compiled program and the time taken in 
running the code. 

The process gg — > gyZ is a leading order one-loop process. It should be both UV as well as IR 
finite. However, individual diagrams may be UV and/or IR divergent. IR divergence is relevant to 
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only light quark cases. All these singularities are encoded in various scalar integrals. To make UV 
and IR finiteness check on our amplitude we have derived all the required scalar integrals (upto box 
scalar integrals with two massive external legs) analytically following t'Hooft and Veltman [|l3|]. 
We regulate the UV divergence of the scalar integrals using dimensional regularization and infrared 
singularities by using small quark mass (the mass regularization). Our results are in agreement with 



those given in [ 14]. Following are the details of various checks made on our amplitude, 



1 . UV Finiteness: 

The ultimate source of the UV divergence is the bubble scalar integral, Bo- For both the 
massive and massless quark contributions, we have verified that the amplitude is UV finite. 
Each pentagon diagram is UV finite by itself as expected from naive power counting. The 
three classes of box diagrams are also separately UV finite. 

2. IR Finiteness: 

The diagrams with massless internal quarks have mass singularities. Even with small quark 
mass (like bottom quark) these diagrams may have large logarithms which should cancel for 
the finiteness of the amplitude. There can be log 2 (m 9 ) and log(m q ) type of singular terms. 
We have checked explicitly that such terms are absent from the amplitude. We have verified 
that the IR finiteness holds for each box and pentagon diagrams Jl5|]. 

3. Gauge Invariance: 

The amplitude has gauge invariance with respect to three gluons, the photon and the Z- 
boson. This has been checked by replacing the polarization vector of any of these gauge 
particles by its momentum (e^(pi) — > pf) which makes the amplitude vanish. As one would 
expect the pentagon and three classes of box contributions are separately gauge invariant with 
respect to the y and Z. For each gluon one of the three classes of box amplitudes is separately 
gauge invariant and further cancellation takes place among pentagon and the other two box 
contributions. 

4. Decoupling of heavy quarks: 

As a consistency check we have also verified that the vector part of the amplitude vanishes 



in the large quark mass limit Jig]. This feature of the amplitude is very closely related to it's 
UV structure. The decoupling holds for each penta amplitude and also for each class of box 
amplitudes. 



4. Numerical Results 



We have computed the amplitude numerically using the real polarization vector basis. How- 
ever, given the number of diagrams, the number of polarization combinations and the length of the 
amplitude, we have to run the code in parallel environment using a PVM implementation of the 
VEGAS algorithm [|l7|, [18[ ]. We have used more than 20 cores to run the code in parallel environ- 
ment. Still it takes more than a day time to get suitable cross-section. To obtain the distributions of 
the kinematic variables, one will have to run the code on more cores for a longer time. 
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Figure 4: Variation of the cross-section (a) with the centre of mass energy 



In such loop calculations, there is an issue of numerical instabilities due to large cancellations 
in the evaluation of the scalar and tensor integrals. To minimize this problem, we have used the 
OneLoop implementation of the scalar integrals, So, Co, and Do. Even with this implementation, 
we have faced numerical instabilities in the evaluation of pentagon-tensor integrals. This is related 
to inaccurate evaluation of Gramm determinants for certain phase-space points. We have used stiff 
cuts on the pj to reduce the impact of this problem. With this, we have about 1 in 10,000 points 
that give rise to numerical instabilities. We need to study this issue more closely. For the time 
being we judiciously throw away such points. This is a general practice. Since such points are few, 
given the nature of the integral evaluation, we don't expect significant impact of this aspect of the 
calculation. 

In Fig (Q), we present the results of the cross-section calculation. These results include fol- 
lowing kinematic cuts: 

Pp ZJ > 50GeV, |Tf' Z ' j | < 2.5, R(y,j) > 0.6. 

We have also chosen factorization and renormalization scales as \ij = = . For practical rea- 
sons, we have used CTEQ4 NLO parton distribution functions. (This is because of the availability 
of the parametrized versions of the parton distribution functions without the external tables.) We 
note that at typical LHC energy, the cross-section would be of the order of 40 fb. Still one may 
expect few thousand of such events after a few years of LHC operation at 14 TeV CM energy. If 
we relax the cut on the pj of the photon and the jet, then the cross-section would increase signif- 
icantly. Our cut is stiff to avoid the numerical instability. As we discussed above the axial-vector 
contribution will also increase the cross-section. 
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5. Conclusions 

We have presented the results of the cross-section calculations for the process gg — >■ gyZ. This 
process occurs at one-loop level via pentagon and box diagrams. We have made a number of 
checks on our calculation. We have verified the cancellation of UV and mass singularities. We 
have also checked gauge invariance with respect to all the gauge particles. We have plotted the 
cross-section as a function of the centre of mass energy. The results are for the vector contribution 
only. This contribution is gauge invariant. Axial-vector contribution does not interfere with the 
vector contribution and it is separately gauge invariant. This axial-vector contribution can only 
increase the cross-section. We have obtained the cross-sections using the parallel running of the 
code. With a stiff cut of about of 50 GeV on the transverse momentum of the photon and the jet, 
the cross-section is about 40 fb at the centre of mass energies of 14 TeV. 
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